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Abstract 

Pearson and Bellissard recently built a spectral triple — the data of Riemanian noncommuta- 
tive geometry — for ultrametric Cantor sets. They derived a family of Laplace-Beltrami like 
operators on those sets. Motivated by the applications to specific examples, we revisit their 
work for the transversals of tiling spaces, which are particular self-similar Cantor sets. We use 
Bratteli diagrams to encode the self-similarity, and Cuntz-Krieger algebras to implement it. 
We show that the abscissa of convergence of the (^-function of the spectral triple gives indica- 
tions on the exponent of complexity of the tiling. We determine completely the spectrum of 
the Laplace-Beltrami operators, give an explicit method of calculation for their eigenvalues, 
compute their Weyl asymptotics, and a Seeley equivalent for their heat kernels. 
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1 Introduction and summary of the results 



In a recent article [36] , Pearson and Bellissard defined a spectral triple — the data of Ricmanian 
noncommutative geometry (NCG) [11] — for ultrametric Cantor sets. They used a construction 
due to Michon [35] : any ultrametric Cantor set (C, d) can be represented isometrically as the set 
of infinite paths on a weighted rooted tree. The tree defines the topology, and the weights encode 
the distance. The spectral triple is then given in terms of combinatorial data on the tree. 

With this spectral triple, they could define several objects, including a ^-function, a measure /i, 
and a one-parameter family of operators on L^{C, ^), which were interpreted as Laplace-Beltrami 
operators. They showed the abscissa of convergence of the (^-function to be a fractal dimension of 
the Cantor set (the upper box dimension). 

The goal of Pearson and Bellissard was to build a spectral triple for the transversals of tiling 
dynamical systems. This opened a new, geometrical approach to the theory of tilings. Until 
now, all the operator-algebraic machinery used to study tilings and tiling spaces was coming from 
noncommutative topology. Striking application of noncommutative topology were the study of 
the iiT-theory for tiling C*-algebras [3,5,30,31], and namely [1,13,15] for computations applied 
to substitution tilings. A follow-up of this study was the gap-labeling theorems for Schroedinger 
operators [4,6,7,25,42] — a problem already appearing in some of the previously cited articles. 
Other problems include cyclic cohomology and index theorems (and applications to the quantum 
Hall effect) [2, 8, 26-29, 32]. .. These problems were tackled using mainly topological techniques. 
The construction of a spectral triple is a proposition for bringing geometry into play. 

In this article, we revisit the construction of Pearson and Bellissard for the transversals of 
some tiling spaces. For this purpose, we use the formalism of Bratteli diagrams instead of Michon 
trees. This approach is equivalent, and applies in general to any ultrametric Cantor set. But in 
some cases, the diagram conveniently encodes the self-similarity, and is very well suited to handle 
explicit computations. For this reason, we will focus on diagrams arising from substitution tilings, 
a class of tilings which is now quite well studied [1,33,38,41]. 

Bratteli diagrams were introduced in the seventies for the classification of Ai^-algebras [9]. 
They were used by Versik to encode measurable Z-actions, as a tool to approach some dynamical 
systems by a sequence of periodic dynamical systems [43] . They were then adapted in the topolog- 
ical setting to encode Z-actions on the Cantor set [13,20,22]. Then, Bratteli diagrams were used 
to represent the orbit equivalence relation arising from an action of Z^ [18] (and recently Z'' [19]) 
on a Cantor set. Yet, it is not well understood how the dynamics itself should be represented on 
the diagram. The case we will look at is when the Cantor set is the transversal of a tiling space, 
and the action is related to the translations. The idea of parametrizing tilings combinatorialy 
dates back to the work of Griinbaum and Shephard in the seventies, and the picture of a Bratteli 
diagram can be found explicitely in the book of Connes [10] for the Penrose tiling. However, it 
took time to generalize these ideas, and to understand the topological and dynamical underlying 
questions. 
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Figure 1: A self-similar Bratteli diagram associated with the matrix (root on the left). 

In the self-similar case (for example when the Cantor set is the transversal of a substitution 
tiling space), the diagram only depends on an adjacency (or Abelianization) matrix. There is a 
natural C*-algebra associated with this matrix, called a Cuntz-Krieger algebra [12]. Its generators 
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implement recursion relations and therefore provide a method of eomputation for the eigenvalues 
of the Laplace-Beltrami operators. 

While Bratteli diagrams are suited to facilitate computations for any self-similar Cantor set, 
we focus on diagrams associated with substitution tilings. Indeed, a tiling space comes with a 
convenient distance, which is encoded (up to Lipschitz-equivalence) in a natural way by weights 
on the diagram. 

Results of the Paper 

Let S be a weighted Bratteli diagram f Definition 12.21 and 12. 9p . and [dB.dw) the ultrametric 
Cantor set of infinite paths in B (Proposition I2.10[) . 

The Dixmier trace associated with the spectral triple gives a measure ^Y>ix on dB. The 
construction of Pearson-Bellissard gives a family of Laplace-Beltrami operators Ag, s G R, on 
L^(9;B, d/iDix), see Definition 14.121 For all s, is a nonpositive, self-adjoint, and unbounded 
operator. For a path 7, we denote by [7] the clopen set of infinite paths with prefix 7, and by x-y 
its characteristic function. And we let exti(7) denote the set of ordered pairs of edges that extend 
7 one generation further. The operator A^ was shown to have pure point spectrum in [36]. In this 
paper, we determine all its eigenelements explicitly. 

Theorem 14.31 The eigenspaces of Ag are given by the subspaces 

= ( 1 TX-r-a 1 — jiX-f-b ■■ {a,b) G exti(7)\ 

for any finite path 7 in B. We have dhnE-y ~ — 1 where is the number of edges in B 
extending 7 one generation further. 

The associated eigenvalues A7 are also calculated explicitly, see equation (|4.13p . An eigenvector 
of As is simply a weighted sum of the characteristic functions of two paths of the same lengths 
that agree apart from their last edge, see Figure [5] for an example. 

o o 




Figure 2: Example of eigenvectors of A3 for the Fibonacci diagram. 

The transversal S to a substitution tiling space (fi, oj) of M.'^, can be described by a stationary 
Bratteli diagram, like the Fibonacci diagram shown in Figure [T] There is a "natural map", called 
the Robinson map in [30], ijj : dB — *■ 2, which under some technical conditions (primitivity, 
recognizability, and border forcing) is a homeomorphism fTheorem l2.22p . We endow 2 with the 
combinatorial metric ds- two tilings are £-close if they agree on a ball of radius 1/e around the 
origin. Let A be the Abelianization matrix of the substitution, and App its Perron-Frobenius 
eigenvalue. We denote by ^(7) the weight of a finite path 7 in S (Definition l2.9|l . 
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Theorem 12.251 // there are constants c+ > c_ > such that c_Ap^' < w{'j) < c^Ap^ for 
all paths 7 of lengths n, then the homeomorphism ip : dB S is (d^-ds.) bi-Lipschitz. 

In our case of substitution tiling spaces, we have a (^-function which is given as in [36] by 

as) = ^ w{^r , 

where IT is the set of finite paths in B. It is proven in [36] that, when it exists, the abscissa of 
convergence so of ^ is the upper box dimension of the Cantor set. For self-similar Cantor sets it 
always exists and is finite. 

Theorem 13.71 For a weighted Bratteli diagram associated with a substitution tiling space of di- 
mension d, the abscissa of convergence of the (^-function is sq ~ d. 

We also have an interpretation of Sq which is not topological. We link sq to the exponent of 
the complexity function. This function p, associated with a tiling, counts the number of distinct 
patches: p(n) is the number of patches of radius n (up to translation). We present two results. 

Theorem 13.121 For the transversal of any minimal aperiodic tiling space with a well-defined 
complexity function, the box dimension, when it exists, is given by the following limit: 

~ 1- lnb(".)) 
dmi ^ = Imi sup — — — — — . 
n-^+oo ln(n) 

And we deduce the following. 

Corollary 13.141 Let S be the transversal of a substitution tiling of dimension d, with complexity 
function p. Then there exists a function v such that: 

p{n) = n'^^'^\ with lim iy{n) = d. 

n — >-|-oo 

With the above choice of weights we can compute the Dixmier trace iiuix- Furthermore, there 
is a uniquely ergodic measure on (r2,M''), which was first described by Solomyak [41]. It restricts 
to a measure fi^^ on the transversal, and we have the following. 

Theorem 13.81 With the above weights one has 4'*ifJ-Dix) = fJ-vf^ ■ 

As B is stationary, the sets of edges between two generations (excluding the root) are isomor- 
phic. Let us denote by £ this set, and by Sq the set of edges linking to the root. Thanks to the 
self-similar structure we can define affine maps u^, e G £, that act on the eigenvalues of as 
follows (see section [5TT|) 

where Uerj is an extension of the path r/ (see Definition l5.16p . and (3^ a constant that only depends 
on e. Those maps Ue,e G £, form a representation of a Cuntz-Krieger algebra associated with 
the matrix A. If 7 = (e, ei, 62, ■ ■ ■ e„), e G Ci G is a path of length n, let be the map 
o o ■ • 'Ue,^. Let US dcnotc by X^s G So-, the eigenvalues of corresponding to paths of 
length 1. For any other eigenvalue A-y, there is a (unique) such that 

n 

X^ = U^iXe) = A^Ae + ^ Ai-l/3e, , 

J = l 

where A^ = A^pp^ ^ That is, the Cuntz-Krieger algebra allows to calculate explicitly the full 
spectrum of A^, from the finite data of the A^, e G £q, and /3e7 e G £ — which are immediate to 
compute, see Section [O] 

For instance, for the Fibonacci diagram (Figure [T]) and s ^ sq ^ d, there are only two such 
maps, namely Ua {x) = x - (jP — 0, and Uh{x) = x ■ (p"^ -\- cj), where ip — {1 + ^/5) /2 is the golden mean. 
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Figure 3: Distribution of the eigenvalues for the Fibonacci diagram. 

The eigenvalues of are all of the form p + q(jr' for integers p,q Z. They can be represented as 
points (p, q) in the plane; these points stay within a bounded distance to the line directed by the 
Perron-Frobenius eigenvector of the Abelianization matrix, see Figure [31 This is an example of a 
general result, valid for hyperbolic substitutions, see Theorem 15.71 We treat further examples in 
Section [6l in particular the Thue-Morse and Penrose tilings. 

The name Laplace-Beltrami operator for can be justified by the following two results. For 
s = sq = d, in analogy with the Laplacian on a compact d-manifold, A^ satisfies the classical 
Weyl asymptotics, and the trace of its heat kernel follows the leading term of the classical Seeley 
expansion. 

Let M{X) be the number of eigenvalues of A^ of modulus less than A. 
Theorem 15.51 There are constants < c_ < c+ such that as X ^ +oo one has 

c^X"'^ <Af{X) < c+X'''^. 
Theorem 15.61 There are constants < c_ < c+ such that as t [ Q one has 

Acknowledgements. The authors were funded by the NSF grants no. DMS-0300398 and no. DMS- 
0600956 of Jean Bellissard, while visiting the Georgia Institute of Technology for the Spring term 
2009. It is a pleasure to thank J. Bellissard for his many supportive, insightful, and enthusiastic 
remarks on this work. The authors would like to thank Ian Putnam for useful discussions (which 
lead in particular to the study of the eigenvalue distribution in section [5. 5p . as well as Johannes 
Kellendonk and Marcy Barge for inviting them to Montana State University, and useful discussions 
they had with them there. J.S. also aknowledges financial support from the SFB 701, Universitat 
Bielefeld, and would like to thank Michael Baake for his generous invitations. 

2 Weighted Bratteli diagrams and substitutions 

We first give general the definitions of Cantor sets, Bratteli diagrams, and substitution tilings. We 
then describe how to associate diagrams to tilings. 
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Definition 2.1. A Cantor set is a compact, Hausdorff, mctrizablc topological space, which is 
totally disconnected and has no isolated points. 

An ultrametric d on a topological space X is a metric which satisfies this strong triangle 
inequality: 

Vx, y,z <E X, d{x, y) < max{d{x, z), d{y, z)}. 

2.1 Bratteli diagrams 

To a substitutive system, one can naturally associate a combinatorial object named a Bratteli 
diagram. These diagrams were first used in the theory of C*-algebras, to classify AF-algebras. 
Then, it was mainly used to encode the dynamics of a minimal action of Z on a Cantor set. 

Definition 2.2. A Bratteli diagram is an oriented graph defined as follows: 

13 = iV,£tot, r,s), 

Where V is the set of vertices, £tot is the set of directed edges, and r, s are functions £tot V 
(range and source), which define adjacency. We have a partition of V and £tot in finite sets: 

V = IJ V„ ; Stot = U £n, 

n>Q n>0 

Vq is a single element called the root and noted o. The edges of £„ have their source in V„ and 
range in V„+i, that is: 

r : £n ^ Vn+i, s : 
We ask that s~^{v) ^ for all veV, and r-^{v) 7^ for aU v e V \ Vq. 

Remark 2.3. If, for all w e V \ {o}, r^^{v) is a single edge, then S is a tree. In that sense, the 
formalism of Bratteli diagrams includes the case of trees, and so is a generalization of the case 
studied in [36]. However, our goal is to restrict to self-similar diagrams, for which computations 
are easier. 

Definition 2.4. A path 7 of length n G N U {+00} in a Bratteli diagram is an element 

n-l 

(£0,61,62, ...) e Y[ 

which satisfies: 

for all < i < n, ''(e,;) ~ s(ei+i). 

We call n„ the set of paths of length < +00, 11 the set of all finite paths, and dB the set of 
infinite paths. 

The function r naturally extends to 11: if 7 = (eg, ■ • ■ , 6„), then r(7) := r(e„). 
In addition to the definition above, we ask that a Bratteli diagram satisfies the following 
condition. 

Hypothesis 2.5. For all v G V, there are at least two distinct infinite paths through v. 

Remark 2.6. If, for all n > 1 and (v,v') S Vn x Vn+i, there is at most one edge from v to w', we 
can simply encode path by vertices: the following map is an homeomorphism onto its image. 

(eo,6i, . . .) I — > (eo,r(6i), . . .). 

Definition 2.7. Given two finite or infinite paths 7 and 7', we note 7 A 7' the (possibly empty) 
longest common prefix of 7 and 7'. 
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The set dB is called the boundary of B. It has a natural topology inherited from the prod- 
uct topology on Jlit^^«' which makes it a compact and totally disconnected set. A basis of 
neighborhoods is given by the following sets: 

[7] = {x G dB ; 7 is a prefix of x}. 

Hypothesis 12.51 is the required condition to make sure that there are no isolated points. This 
implies the following. 

Proposition 2.8. With this topology, dB is a Cantor set. 

Definition 2.9. A weight on dB is a function w '.V ^ M^, satisfying the following conditions: 

(i) w{o) = 1; 

(ii) sup{z/;(v) ; v € V„} tends to when n tends to infinity; 

(iii) Ve G ftot, w{s{e)) > w{r{e)). 

A weight extends naturally on paths: by definition, w(7) := w{r('y)). 
Proposition 2.10. We define a function on (dB)^ by: 

w{r{xAy)) if x ^ y; 
otherwise 

It is a ultrametric on dB, which is compatible with the topology defined above, and so {dB,dw) is 
a ultrametric Cantor set. 

A case of interest is when the Bratteli diagram is self-similar. For a diagram, self-similarity 
means that all the V„ are isomorphic, and all £„ are isomorphic for n > 1, and r, s commute 
with these identifications. We will focus on self-similarity when the diagram is associated with a 
substitution. 

2.2 Substitution tilings 

Let us give the definition of a tiling of . The tilings we are interested in are constructed from 
a prototile set and an inflation and substitution rule on the prototiles. The notion of substitution 
dates back to the sixties. For example the self-similarity of the Penrose tiling already appeared 
in [37]. However, systematic formalisation of substitution tilings and properties of their associated 
tiling space was done in [33, 40] in the nineties. In this section, we follow the description by 
Anderson and Putnam [1] with some minor changes. In particular, some non trivial facts are cited 
along the text. We do not claim to cite the original authors for all of these. The reader can refer 
to the reviews [39] and [16]. 

By tile, we mean a compact subset of W^, homeomorphie to the closed unit ball. A tile is 
punctured if it has a distinguished point in its interior. A prototile is the equivalence class under 
translation of a tile. Let ,4 be a given set of punctured prototiles. All tilings will be made from 
these tiles. 

A patch p on A\s a. finite set of tiles which have disjoint interior. We call A* the set of patches 
modulo translation. A partial tiling is an infinite set of tiles with disjoint interior, and the union 
of which is connected, and a tiling T is a partial tiling which covers , which means: 

Supp(r) := IJ t = K'^ . 

teT 

A substitution rule is a map w which maps tiles to patches, and such that for all tile p €z A, 
Supp(a-'(p)) = ASupp(p) for some A > 1. The factor A inflates the tiles p, which is then cut into 
pieces; these pieces are elements of A. The map to extends to patches, partial tilings and tilings. 

Since we will represent tilings spaces by diagrams, we have a specific interest for the combina- 
torics of the substitution. 



dw{x,y) 
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Definition 2.11. Given a substitution to, its Abelianization matrix is an integer- valued matrix 
= iapq)p,qeA (or simply A) defined by: 

Vp, q Cz A, Upq = number of distinct translates of p included in uj(q). 

We now define the tiling space associated with uj. 

Definition 2.12. The tiling space is the set of all tilings T such that for all patch p of T, p is 
also a subpatch of some a;"(i) {n E N, t E A). 

We make the following assumptions on uj: 

Hypothesis 2.13. (i) The matrix A is primitive; A"' has non-negative entries for some n. 

(a) Finite local complexity (FLC): for all R > 0, the set of all patches p d T which can be 
included in a ball of radius R is finite up to translation, for all T E O. 

Just as every tile has a distinguished point inside it. Q has a distinguished subset. This is the 
Cantor set which will be associated with a Bratteli diagram. 

Definition 2.14. Let S be the subset of il of all tilings T such that is the puncture of one of 
the tiles of T. This set is called the canonical transversal of fi. 

Note that, by this definition, the substitution lo extends to a map ^ fi. 

We will assume in the following that w : ^ is one-to-one. This condition is equivalent 
(Solomyak [40]) to the fact that no tiling in has any period. Furthermore, it implies that is 
not empty, and w : f2 — > is onto. 

The set of interest for us is S. Its topology is given as follows. For any patch p, define the 
following subset of S: 

Up = {TEE; pcT}. 

Note that Up can be empty. Nevertheless, the family of all the C/p's is a basis for a topology on S. 
Let us now define a distance d on Q: 

d{T, T') = inf (^{e > ; 3a;, y E B{0, e) such that 

i?(0, 1/e) C Supp ((T -x)n {T' -y))}u {1}^ . 

Two tilings are d-close when, up to a small translation, they agree on a large ball around the 
origin. This distance, when restricted to S, is compatible with the topology defined above. With 
this topology, S is a Cantor set. 

Furthermore, the map a; : f2 ^ is a homeomorphism, and the dynamical system (r2,R'') 
given by translation is continuous and uniquely ergodic [41]. 

Substitution tiling spaces are minimal, which means that every R''-orbit is dense in fl. Com- 
binatorially, this is equivalent to the fact that these tilings are repetitive: for all R> 0, there is a 
bound PR, such that every patch of size R appears in the tiling whithin range pR to any tile. 

We need an additional assumption — the border forcing condition. It is required in order to 
give a good representation of S as the boundary of a Bratteli diagram. 

Definition 2.15. Let p E Che a patch. The maximal unambiguous extention ofp is the following 
patch, obtained as an intersection of tilings: 

Ext(p) = fl {t e f7 ; per}. 

So any tiling which contains p also contains Ext(p). The patch Ext(p) is called the empire of 
the patch p by some authors [21]. 
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Definition 2.16. Assume there is some m e N and some e > such that for all t £ A, 

Ext contains an el-ncighborhood of uj''^{t). 

Then lo is said to force its border. 

It is always possible to give labels to the tiles, in order to change a non-border forcing substi- 
tution into a border-forcing one. See [1]. 

Example 2.17. Let us give an example in dimension one. Consider the substitution defined on 
symbols by: 

a 1-^ baa 
& i~> 6a 

It has a geometric realization, with a being associated with an interval of length (f) = {1 + a/5)/2, 
and b with an interval of length 1. Then the map above is a substitution in the sense of our 
definitions, with X = (p. The transversal S is naturally identified with a subset of {a, 6}^ (the 
subshift associated with to). The map lo satisfies the border forcing condition. Indeed, if p is a 
patch of a tiling T (which we identify symbolically with a finite word on {a, 6} in a bi-infinitc 
word), let x, y G {a, 6} be the letters preceding and following p in T respectively. Then u}{x) ends 
by an a and uj{y) begins by a b. Therefore, w(p) is always followed by b and preceded by a. This 
proves that uj forces its border. 

2.3 Bratteli diagrams associated with substitution tilings 

We show how to associate a Bratteli diagram to a substitution, and identify the transversal with 
the boundary of the diagram. It is clear that these two are homeomorphic, being Cantor sets. We 
will give an explicit and somehow natural homeomorphism. 

Let cj be a primitive, FLC and border forcing substitution, and A = {aij)ij^A its Abelian- 
ization matrix. Let A be the expansion factor associated with lo. Let S be the transversal of the 
tiling space associated with uj. 

The diagram B associated with the substitution is defined as follows. 

Definition 2.18. The diagram associated with uj is the diagram B = {V,£tot,r, s), where: 

Vn > 1, Vn^Ax {n}, 

for all n > 1, for all i,jEA, there are exactly Uij edges of from (i, ji) to {j,n + 1), and for all 
V Cz Vi there is exactly one edge e„ from the root o to v. 

We call £ the set of "models" of edges from one generation to another. All £"„ are copies of £ 
(for example with identification £„ = £ x {n}). Since the models for the vertices are the elements 
of A, we do not use a specific name. The maps source and range can be restricted as maps £ ^ A. 
When the "depth" of a vertex is not important, we will sometimes consider r, s as functions valued 
in A. 

Remark 2.19. Combinatorially, the diagram only depends on the Abelianization matrix of uo. It is 
indeed possible to associate a diagram to a primitive matrix A with integer coefficients. We would 
have similar definitions, with V„ = I x {n}, with / the index set of the matrix. 

Proposition 2.20. This is a Bratteli diagram in the sense of Definition [27Sl except in the case 
A={1). 

Example 2.21. Figure [1] is an example of the self-similar diagram associated with the Fibonacci 
substitution: a t-^ ab; b t-^ a. 

There is a correspondence between the paths on B and the transversal S. It depends on a 
choice on the edges which remembers the geometry of the substitution: each e e £ from a G ^ to 
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b € A corresponds to a different occurrence of a in uj{b). This correspondence is a homcomorphism 
(/) : S ^ dB, called the Robinson map, as defined in [30]. We first give the definition of (f>, and we 
will then give a condition the weight function w so that is bi-Lipschitz. 

To construct (f), start with T G S. Let t e T be the tile containing 0, and [t] G A the 
corresponding prototile. Then, the first edge of is the edge from o to ([i], 1) G Vi. Assume 
that the prefix of length n of (j>{T) is already constructed and ends at vertex where t' 

is the tile of (jj~'^"~-'^'(T) which contains 0. Let t" be the tile containing in (jj~"(r); then the 
(n + l)-th edge of (^(T) is the edge corresponding to the inclusion of t' in Lo{t"). This edge ends 
at ([t"], n + 1). By induction we construct 4>{T). 

Theorem 2.22 (Theorem 4 in [30]). The function (/) :'E. ^ dB is an homcomorphism. 

We give explicitly the inverse for as follows. Let us define : H — > A* which is increasing 
in the sense that if a path is a prefix of another, then the patch associated to the first is included 
in the patch associated to the second. The image of an infinite path will then be defined as the 
union of the images of its prefixes. 

If (e) is a path of length one, then -0(e) is defined as the tile r{e) with puncture at the origin. 
Now, given a path 7 of length n, assume its image by -0 is well defined, and is some translate 
of a'"~^(r(7)), with the origin at the puncture of one of its tiles. Consider the path 7.e, with 
e G £„ such that s(e) — ^(7). Then, e encodes an inclusion of s(e) inside uj{r{e)). It means that it 
encodes an inclusion of w"~^(r(7)) inside [j"(r(e)). The patch ipil-^) is defined as the translate 
of [j"(r(e)) such that the inclusion ^(7) C tpij-e) is the inclusion defined by the edge e. 

Remark that if x G dB, ipix) does not a priori define more than a partial tiling. The fact 
that it corresponds to a unique tiling results from the border forcing condition. Then the map 0, 
extended to infinite paths, is the inverse of (f>. 

Continuity can be proved directly, but can also be seen as a consequence of Theorem I2.25| 
which we prove later. 

Lemma 2.23. There exists Ci, C2 positive constants, such that for all path 7 of length n in B, 



where A is the expansion factor ofuj. 

Proof. Let k be the smallest number such that for all w G V, there exists two distinct paths of 
length k starting from v. For example, k = 1 when there are two elements of Stot starting from 
V, for all V, and fc = 2 in the case of Fibonacci, pictured in Figure [TJ Let C be the maximum 
diameter of the tiles. Let 7 be a path of length n. Then, one can find two distinct extensions of 
length n + fc of 7; call them 71 and 72. Remember that ^(7^) is some translate of u!"^^~^(ti) with 
ti G A, and G Supp (0(7)). So ■!/'(7i) ^-nd ^'(72) differ within range CA"+'^^^. Therefore, with 
Ci = CA*^"^, we have the second line. 

For the first inclusion, let m be the exponent for which uj satisfies the border forcing condition, 
and 7 a path of length m. Let i G -4 be the range of 7. Then Ext{uj™' (t)) covers an e-neighborhood 
of Supp(ijj"'(t)). It means that, since "0(7) is a translate of a;™(i) which contains 0, 



With n = m + k and Ci = e/A™, one has the result. For n < m, the inequality still holds (up to 




i?(0,e)cExt(^(7)). 



Similarly, for all path of length m + k 



B{0,e\'') C Ext(^(7)). 



a reduction of e). 



□ 
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Remark 2.24. Wc can simplify the diagram in the ease where a symmetry group G C 0„(M'') acts 
freely on A, such that 

Vg e G, Vp e A, uj{g ■ p) ^ g ■ uj{p). 

With this, it is possible to extend the action of G to the edges of £tot- It induces naturally an 
action of G on dB. The group G also acts on 51 by isometrics, and these two actions are conjugate 
by ip. 

One can "fold" the diagram by taking a quotient as follows. Define: B' — {V' ,£[a^,r' , s'), 
where Vq = V, = £o x G, and 

V?! > 1, V; = Vn/G and £'„ = f„/G. 

Furthermore, for [ejc e £n/G, r'{[e]G) := [''(e)]G and s'([e]G) [s(e)]G- These definitions do not 
depend on the choice of the representant e. For (e,g) S £q = x G, define r'{{e,g)) := [r{e)]G 
and s'{{e,g)) := o. One can check that these two diagrams are "the same" in the sense that the 
tree structure of their respective sets of paths of finite length is the same. 

In terms of substitution, the image by of a path 7 = (eq, ei, . . . , e„_i) is the image of 
a)"(r(7)) under some element of R'' xi G . The translation part is given by the truncated path 
(ei, . . . , e„_i). The rotation part is encoded in the first edge Eq G fg- ^0 ~ (e', 9), where e' brings 
no additional information, but g corresponds to the choice of an orientation for the patch. 

There is still an action of G on dB', defined by: g ■ ((eq, h),ei, . . .) = {{eo,gh), ei, . . .) . 



2.4 Weights and metric 

We gave in Theorem 12.221 an explicit homeomorphism (/) between the boundary of the Bratteli 
diagram, and the transversal of the tiling space. Since we are interested in metric properties of the 
Cantor set, we now show that a correct choice of weights on the vertices of the Bratteli diagram 
gives a metric on dB which is Lipschitz equivalent to the usual metric on S (the equivalence being 
induced by </>). 

Theorem 2.25. Let B the Bratteli diagram associated with a substitution lo . Let A be the inflation 
factor of to. We make the following assumption on the weight function w: 

Vn > 1, Vw € V, w{v, n + 1) = —w{v, n). (2-1) 

A 

Then the function (j) : (S, d) — > {dB,d.^) defined in Provosition \2.2^ is a bi-Lipschitz homeomor- 
phism. 

Proof. It is enough to show that there are two constants m, M > 0, such that for all 7 € 11, 

m < '^^^"^(^^^ < M. 
- diam((/)-i([/^)) - 

Since the are a basis for the topology of dB, this will prove the result. 

Let 7 G n„. By Lemma [2.231 any two tilings in (p~^{U^) coincide on a ball of radius at least 
GiA". Therefore, 

diam((/.-i([/^)) < -^A-". 

On the other hand, it is possible to find two tilings in (j)'~^(U^), which disagree on a ball of radius 
G2A". Therefore, 

diam((/)-i(t/.^)) > T^A-". 

And by definition of the weights, 

min{u;(w, 1) ; w £ V}A-"+i < diam(C/^) < max{u)(u, 1) ; w G V}A-"+^ 
Together with the previous two inequalities, this proves the result. □ 
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3 Spectral triple, (^-function, and complexity 



3.1 Spectral triple 



Let S be a weighted Bratteh diagram, and let {dB, d) be the uhrametric Cantor set of infinite 
rooted paths in B. Pearson and BeUissard built in [36] a spectral triple for {dB, d) when S is a tree 
(that is when for all vertex v G V \ {o}, the fiber r~^(v) contains a single point). In our setting, 
their construction is adapted as follows. 
A choice function on S is a map 



n 

7 



dBxdB 



such that (i(r+(7), t_(7)) ~ diam[7] , 



(3.2) 



and we denote by E the set of choice functions on B. Let Cup{dB) be the C*-algcbra of Lipschitz 
continuous functions on {dB, d). Given a choice t € E we define a faithfuU ^-representation tt^ of 
Cup{dB) by bounded operators on the Hilbert space H = P{H) 



as 



■yen 



f{r+il)) 
/(t-(7)) 



(3.3) 



This notation means that for all ^ € 7i and all 7 G LI, 

(^r(/)-e)(7) 

A Dirac operator D on Ti. is given by 







■e(7). 



^ diam[7] 



1 

1 



that is D is a self-adjoint unbounded operator such that {D^ 
tator 

/(T+(7))-./(r_(7)) 



(3.4) 



1) ^ is compact, and the commu- 



[A -.(/)] =0 
7en 



diam[7] 



-1 

1 



(3.5) 



is bounded for all / G C'up (dB) . Finally a grading operator is given by F = 1 



;2(n) 



I 
-1 



and satisfies F^ = F* = F, and commutes with tTt and anticommutes with D. The following is 
Proposition 8 in [36]. 

Proposition 3.1. (CLip(<9S), ?i, tTt-, D, F) is an even spectral triple for all t E E. 

In [36] the set of choice functions E is considered an analogue of a tangent bundle over dB, 
so that the above commutator is interpreted as the directional derivative of / along the choice 
T. The metric on dB is then recovered from the spectral triple by using Connes formula, i.e. by 
taking the supremum over all directional derivatives 



Theorem 3.2 (Thm. 1 in [36]). The following holds: 

d{x,y)=snp{\f{x)-f{y)\ ; f e Cup{dB) , sup || [i?, 7r,(/)]|| < l} 

reE 

Definition 3.3. The (-function associated with the spectral triple is given by: 

1. 



C(s) 



-Ti- \D\ 



yen 



(3.6) 



(3.7) 



and we will denote by sq G M its abscissa of convergence, when it exists. 
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We now assume that the weight system on B is sueh that sq G M. 

Definition 3.4. The Dixmier trace of a function / e Cup{dB) is given by the following limit, 
when it exists: 

^^■'^ slso TrdDl-^) ^ ^ 

It defines a probability measure on dB and does not depend on the choice t E E (Theorem 3 
in [36]). Furthermore, when / = X7 is a characteristic function, the limit above can be rewriten: 



hm ^ 

where 11^ stands for the set of paths in 11 with prefix 7 



3.2 The measure on dB 

Let (fi,cj) be a substitution tiling space, A be the Abelianization matrix of cj, and B be the 
associated Bratteli diagram, as in previous section (Definition I2.18P . We assume that B comes 
together with a weight function w, which satisfies the properties given in Definition 12.91 and 
adapted to the substitution uj. In particular, it satisfies w(a,n + 1) = \^^w{a,n) for all a E A 
and n € N. 

As we assumed A to be primitive, it has a so called Perron-Frobenius eigenvalue, denoted 
Api?, which satisfies the following (see for example [23]): 

(i) hpF is strictly greater than 0, and equals the spectral radius of A; 

(ii) For all other eigenvalue v of A, \v\ < App; 

(iii) The right and left eigenvectors, vp and v^, have strictly positive coordinates; 

(iv) If vp and vl are normalized so that {vp, vp) — 1, then: 

A" 

lim = wlVr; 



(v) Any eigenvector of A with non- negative coordinates corresponds to the eigenvalue App. 

A classical result about linear dynamical systems together with the properties above gives the 
following result, which will be needed later. 

Lemma 3.5. Let A be a primitive matrix with Perron-Frobenius eigenvalue App. Let Pm be its 
minimal polynomial: Pm ~ {X — App) J^f^j^ {X — a^)™^*''. Then, the coefficients of A"- are given 
by: 

[A'%b = CabAlp + ^/"'VX , 
i=l 

where the Pi 's are polynomials of degree m[i), and Cab > 0. 

Note that the coefficients (a, b) of A" gives the number of paths of length n in the diagram 
between some vertex (a, k) G Vk and (&, n + k). This lemma states that this number is equivalent 
to CabApp when n is large. 

Proof. We have Cab > 0, as it is the (a, b) entry of the matrix vpv]^ defined above. The rest is 
classical, and results from the Jordan decomposition of the matrix A. □ 

We assume that $7 is a c?-dimensional tiling space. Since co expands the distances by a factor 
A, the volumes of the tiles are dilated by A'^. This gives the following result: 
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Proposition 3.6. Let App be the Perron-Frohenius eigenvalue of A. Then App = X'^ . In 
particular, App > 1. 

Theorem 3.7. The (^-function for the weighted Bratteli diagram [dB,w) has abscissa of conver- 
gence sq — d. 

Proof. We have 

c(s) = ^«;(7r- 

The quantity w{'^) only depends on r(7). Furthermore, w{a,n) tends to zero hke A" = (Ap/^)" 
when n tends to infinity. So we have: 

J2 ™"Ap/''Card(n„) < C(s) < ^ A/^Ap/'' Card(n„), 

neN) neN) 

where m (resp. M) is the minimum (resp. the maximum) of the w{a, 1), a € ^ Now, since Card(n„) 
grows like App up to a constant (sec Lemma 1575]) . wc have the result. □ 

Theorem 3.8. The measure fi given by the Dixmier trace (Definition p.4p ) is well defined, 
and given as follows. Let v = {va)aeA be the (right) eigenvector for A, normalized such that 
See£o " ^- ^'^''^ '^^^ 7 G n, let (a, n) := r{"f) e V„. Then: 

M([7]) = ivAp^+^ 

In particular, in the case of a substitution tiling, ip*{lJ-) the measure given by the frequences 
of the patches, and therefore is the restriction to S of the unique ergodic measure on {Q,R'^) 
(see [41]). 

Proof. Let 7 = {eo,ei, . . . ,e„_i) G n„, and (a,n) :— r(7). Let YL^ the subset of n of all paths 
which have 7 as a prefix. Define: 

J., , »?en^ 

Then /^([7]) — limsjd/(s), when this limit exists. The terms of the sum above can be grouped 
together: if 77 is a path of H^, then the quantity w{riy only depends on the length of 77, say n + k 
{k > 0), and on r(ri) = b G A. Then, if we call N{a, b; k) the number of paths of length k from a 
to 6, we can group the sum and write: 

^^7V(a,6;/c)w(6,n + fc)'* 

, fe>OhG.A 



1 + EE E^w^)'^;^M^'^)' 

*:>0 £££0 beA 



Now, since w{b, n)" = A'- 1)^, we can write: 

N{a, b; k)w{b, n + kf = eJ A^W{s), 

where Ea is the vector {6a{b))beA, and W is the continuous vector-valued function s ^—t {w{(b, l)y)beA- 
Similarly, 

Y E k)w{b, ky = A(-'^+i'-^i?T A^VK(s), 

ee£o b£A 
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where Eg is the sum over e G £o of aU E^(^^y 
Now, by lemma 1531 we have: 



i=l 

where the Pi's are polynomial in n for s fixed, with Kpp > \ai\ for all z, and Ca(s) > 0. Further- 
more, the Pi are continuous in s. Similarly, 

p 

i=l 

and since Eg^ is a linear combination of the Ea&, 

C£^Y '^'■(e)- (^■^'^) 

Then, we write: 



p / \ \ nt; /v:--^U Pc^O i — 1 



+ E ce{App/XT + E E ^)(«'VA^)^ 

fc>0 A;>0 i=l 

Eca(Api./A^)" + i?i(s) 



fc>0 



A- + Ec£(ApF/A«r + i?2(s) 

fc>0 

Note that the expression above is defined a priori for s > d, but i?i {i = 1,2) is defined and 
continuous for s > d (the continuity results from the absolute convergence of the sum). The 
remaining sums above can be computed explicitly, and we have: 

Cais)+(l-{App/XT)Rlis) 



c£ + (l - (App/A-^)«) (i?2(s) + A-^) ' 

Then it is now clear that this expression is continuous when s tends to d, and has limit Ap^Ca (d) / C£ (d) . 

Let Ua be defined as Ca/cs for all a. Let us show that u ~ v. First, show that {ua)aeA (or 
equivalently, {ca)aeA) is an eigenvector of A associated with App. We have: 

EgA^Wjd) 

Ca{d) = hm — 

= EaLW{d), 

where L ~ xy^ , with x (resp. y) an eigenvector of A (resp. of A^) associated with App, and 
{x,y) = 1. So Ca is the a-coordinate of LW{d) = {y,W{d)) x, and so is an eigenvector of A 
associated to App. Equation (|3.10p now proves that u has the good normalization, and that 
u = V. □ 

3.3 Complexity and box counting dimension 

Definition 3.9. Let (X, d) be a compact metric space. Then, the box counting dimension is 
defined as the following limit, when it exists: 

dim(X,d) = lim-i^^^, 

where Nt is the minimal number of balls of radius t needed to cover X. 
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Theorem 3.10 (Thm. 2 in [36]). Given an ultrametric Cantor set and its associated (^-function, 
let So be the abscissa of convergence of (. Then: 

So = diia{X, d) . 

This dimension can be linked to complexity for aperiodic repetitive tilings. 
Definition 3.11. The complexity of a tiling T is a function p defined as follows: 

p{n) = Card |g C (T - x) ; x e R'^, B{0, n) C Supp(g) and Vq' C q, B{0, n) (f_ Supp(q')}- 

In other words, p(ri) is the number of patches 5 of T of size n. 

Note that when T is repetitive (for example when T is a substitution tiling) the complexity is 
the same for all the tilings which are in the same tiling space. 

Theorem 3.12. Let D, be a minimal tiling space, and S its canonical transversal, endowed with 
the metric defined in Section Let p be the associated complexity function. Then, for this 

metric, the box dimension is given by: 

dMn,d)^ hm 

ln(n) 

Proof. Let Nt be the number of balls of diameter smaller than t needed to cover S. Let us first 
prove that for all n e N, 

p{n) = Ni/n. 

Let C{n) be the set of all patches of size n, so that p{n) — Card(£(n)). Then, since all the tilings 
of S have some patch of size n at the origin, the set: 



[U, ; qeC{n)] 



is a cover of S by sets of diameter smaller than 1 /n. So p{n) > A^i/„. To prove the equality, assume 
we have a covering of S by open sets {Vi ; i £ /}, with Card(/) < p{n), and diam(T/i) < 1/n 
for all n. Then, in every Vi, we can find some set of the form Uq, q a patch. This allows us to 
associate some patch q{i) to all i G /. We claim that for all such q{i), -6(0, n) C Ext(q(i)), where 
Ext(g) is the maximal unambiguous extension of p, as defined in 12.151 Indeed, if it were not the 
case, diam(J7g) would be smaller than n. Therefore, by restriction, to each i € I, we can associate 
a patch q'{i) of size n, such that J7g'(i) C Vi. Since the {Vi}i^i cover 5, all patches of size n are 
obtained this way. and p{n) < Card(/). Since this holds for all cover, p{n) < Ni/^- 
Now, Nf is of course an increasing function of t. Therefore, 

Nt-i<p{[l/t]) <Nt, 

and so: 

Nt-i ^ p{[l/t]) ^ Nt 



ln(t) - ln([l/t]) - ln(i- !)• 
Letting t tend to zero proves the theorem. □ 

Corollary 3.13. Let S be the transversal of a minimal aperiodic tiling space, with a complexity 
function which satisfies Cin" < p{n) < C2n" for Ci, C2, a > 0. Then: 

dim(S, d) = a. 

Let now consider a tiling space associated to a substitution uj. Let B be the weighted Bratteli 
diagram associated with it. We proved in Section [X^ that for a Bratteli diagram associated with 
a substitution tiling of dimension d, the abscissa of convergence is exactly d. It is furthermore 
true that the box dimension of the transversal S with the usual metric is d; it results from the 
invariance of the box dimension under bi-Lipschitz equivalence, which is proved in [14, Ch. 2.1]. 

Therefore, we can deduce the following result: 
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Corollary 3.14. Let CI be a substitution tiling space satisfying our conditions. Then there exists 
a function v such that: 

p{n) = n^^^\ with lim v{n) ~ d. 

n — >+oo 

Equivalently, for all e > 0, there exists Ci, 6*2 > such that for all n large enough, 

Cin'^-" < p(n) < C2ri'^+' . 

This result is actually weaker than what we can actually expect: in fact, there exist Ci, C2 > 
such that 

Cin'* < p{n) < C2n'^. 

The upper bound was proved by Hansen and Robinson for self-afFine tilings (see [39]). The lower 
bound can be proved by direct analysis for substitution tilings. It would also result from the 
conjecture that any d-dimensional tiling with low complexity (which means p{n)/nf' tends to zero) 
has at least one period (see [34]). 

However, it is still interesting to see how the apparently abstract fact that the abscissa of 
convergence sq equals the dimension gives in fact a result on complexity. 



4 Laplace— Beltrami operator 

Let S be a weighted Bratteli diagram. The Dixmier trace (|3.8p induces a probability measure v 
on the set E of choice functions (see [36] section 7.2.). The following is [36] Theorem 4. 

Proposition 4.1. For all s eS. the bilinear form on L^{dB,dfi) given by 

Qsif,9) = lJ^^H\Dr[D.Mf)r[D.M9)]) dv{T), (4.11) 

with dense domain DomQs = (X7 ■ 7 G n) , is a closable Dirichlet form. 

The classical theory of Dirichlet forms [17] allows to identify Qs{f,g) with (/, A^g) for a 
non-positive definite self-adjoint operator As on L^{dB,dfi) which is the generator of a Markov 
semi-group. We have DomQs C DomA., C DomQs where Qg is the smallest closed extension of 
Qs- The following is taken from [36] section 8.3. 

Theorem 4.2. The operator As is self-adjoint and has pure point spectrum. 

Following Pearson and Bellissard we can calculate As explicitly on characteristic functions of 
cylinders. For a path 77 G H let us denote by exti(77) the set of ordered pairs of distinct edges 
(e, e') which can extend 77 one generation further. 



AsX7 ^ - fr-^((M[7fc] - Ai[7fe+i])X7 - M[7](X7fc - X7fc+i)) (4.12a) 
k=o '^wfcJ 

with Gs (?y) = ^ diam[?7]2-'' ^ ^[77 • e]/i[r; • e'] (4.12b) 

(e,e')Gcxti (r;) 

Note that the term ~ X-yk+i is the characteristic function of all the paths which coincide 
with 7 up to generation k and differ afterwards, i.e. all paths which split from 7 at generation k. 
And fi[jk] — f-ijk+i] is the measure of this set. 

We now state the main theorem which gives explicitly the full spectrum of As. 

Theorem 4.3. The spectrum of As is given by the following, 
(i) is a single eigenvalue with eigenspace (1 = Xob) ■ 
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(a) Aq ~ ITJo) eigenvalue with eigenspace Eq = ij^Xs ^ TTpJ^e' ■ S ^o,^ 7^ £') of 

dimension dinii?o = no — 1, where uq is the cardinality of £q. 



(Hi) For 7 G n, 



A. = 'e (4.13) 



k=0 



is eigenvalue with eigenspace 



= (-r— rX7-e r~~Tr^7-e' '■ (e^e') € exti(7)) (4.14) 

o/ dimension dimii'^ = — 1, where is the number of edges extending 7 one generation 
further. 

Proof. The formula for the eigenvalues is calculated easily noticing that AgX-y-e and A^x-y-e' only 
differ by the last term in the sum in equation (j4.12ap . 

The spectrum of is always the closure of its set of eigenvalues (whatever its domain may 
be). Hence we do not miss any of it by restricting to characteristic functions. 

We now show that all the eigenvalues of are exactly given by the . It suffices to check that 
the restriction of Ag to n„ has exactly dimn„ eigenvalues (counting multiplicity). Notice that x-y 
is the sum of x-y-e over all edges e extending 7 one generation further. Hence an eigenfunction in 

for 7 e Hk can be written as a linear combination of characteristic functions of paths in H„i , for 
any m > k. The number of eigenvalues A-y for 7 S H„_i is '^^^n ^ dim E^ ~ X^yen 1 (""7 ~ 1) ~ 
diniH,! — dimH„_i. So the number of eigenvalues A-y for 7 G H^ for all 1 < fc < n — 1 is 
(diniHfc+i — dimHfc) = dimn„ — dim Hi. And counting and Aq adds up (dim Hi — 1) + 1 
to make the count match. □ 



Remark 4.4. As noted in Remark \2.3\ our formalism with Bratteli diagrams includes as a special 
case the approach of Pearson-Bellissard for weighted Cantorian Michon trees. Hence Theorem \4-.3\ 
gives also the spectrum and eigenvectors of their Laplace- Beltrami operators. 

The eigenvectors (j4.14p are very simple to picture. Given a path 7, and two extensions (a, b) G 
exti(7), an eigenvector is the difference of their characteristic functions weighted by their measures. 
See Figure [5] in section [1] for an example for the Fibonacci diagram. 



5 Cuntz— Krieger algebras and applications 

We now consider stationary Bratteli diagrams. We use the self-similar structure to further char- 
acterize the operator A^, and its spectrum. 



5.1 Cuntz— Krieger algebras 

Let S be a stationary Bratteli diagram. Let A be its Abelianization matrix. Let us denote by £q 
the set of edges e linking the root to generation 1, and by £ the set of edges linking two generations 
(excluding the root). Let A = (ae/)^ be the square matrix with entries Oef = 1 if e can be 
composed with (or followed down the diagram by) / and hef = else. There is an associated 
Bratteli diagram B with Abelianization matrix A, which is "dual" to B in the sens that its vertices 
corresponds to the edges of B and its edges to the adjacencies of edges in B. Note that, because 
the entries of A are zeros or ones, all the edges in B are simple. 

Remark 5.1. If B has only simple edges, we can simply take B ^ B and A = A. 

The Cuntz-Krieger algebra O^, is the C*-algebra generated by the partial isometrics Ue,e € £ 
(on a separable, complex, and infinite dimensional Hilbert space Ti.) that satisfy the following 
relations. 
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^ C*{ Ue, U:,e€£ I UeU:,U:Ue e Vin) , U:Ue = ^f^} ) ' (5.15) 

where V{H) denotes the set of projections in p G V{H) <;=^ = p* = p. 

By abuse of notation we write the basis elements of P(n\ni) as 7 G n\ni. The Cuntz-Krieger 
algebra is represented on ^^(11 \ Hi) as follows. 

jj , X f (e', e, ei, 62, • • • ) if Aeei = 1 Ir \ 

C/e(e,ei,e2,---) = I q ^^^ic (5.16a) 

TT*( \ / (£'.62,63, ■■■) ifei=e /t:r1Gu^ 

[/Je,ei,e2,- • •) = I q gjg^ (5.16b) 

where e' G £0 stands for the (possibly) new edges linking to the root, as illustrated in the case 
of the Penrose substitution below. The orientation of e' however is taken to be the same as 
that of e: if e = (s(ei),.g) then we have e' = (s(e),g) (see Remark I2.24p . In other words we 
require that ipi^' 1 e, ei, • • • ) = w o '^(e, ei, • • • ) + x for some x G M'', and where : OS ^ S is the 
homcomorphism of Proposition 12. 221 See Figure |4] for some examples. 



• • • • • • 




• • • • • • 




Figure 4: Example of some Cuntz-Krieger operators acting on finite paths. 
There is an induced "action" on DomAs C L^{dB,dfi) defined as follows: 

" ^ \ else 

for paths in Iln>2 and by linearity for shorter paths, using the relations x-y = X^e e i^^^ sum 
running over all edges e extending the path one generation further). 

li (fi G Ej, for 7 G n„>2, is an eigenfunction of Ag, then we see from Equation (|5.17p that 
Ue^p, if not zero, is another eigenfunction. Since the diagram B is stationary, we have the following 
obvious fact. 

_ r Eu^^ if C/e7 ^ 



- ^ {0} else: ' " (^-l^) 
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This allows to define an induced "action" of on the set of eigenvalues of Ag as follows. 
Those maps are calculated as shown below. 

Lemma 5.2. Let App be the Perron-Frobenius eigenvalue of A. Let ^ = (e, ei, 62, • • • ) G n„>2. 

Then, if U^J 7^ 0, we have 



Gsio) )^ G,(o) + G«(e') 



Ue(A^)=A},j, ?777^)+ — FTTT^ — + TTTZTx (5-20) 



Proof. By equation (|4.13p we have 



\U 7| — 1 

The terms for = 0, 1, in the above sum give the last two terms in equation (|5.20p . For all 

fc > 2, Ai[(f/e7)fc = Ap],/z[7fe_i], and Gs{{Uel)k) = K^pp''^'''~^Gs{lk-i) (see Theorem EH for the 
rescalling of the measures). Hence the rest of the sum, over k — 2, - ■ ■ \Ue^\ — 1 (and the last term) 
in equation (|5.2ip . rescale by a factor Ap^*^^^^'' to the sum over fc = 1, • • • I7I — 1 (and the last 
term) for the eigenvalue A-y (equation (|4.13p ). We then add the contribution of the root, i.e. the 
term for fc = 0, to get equation (|5.20p . □ 

Note that the maps Ue are affine, with constant terms, written /3e, that only depend on e. We 
will write for now on 

UeiX^) = AA-f + Pe , with A, = A^^+'--''/^ (5.22) 

The eigenelements of corresponding to Eq, and E^, e €z £q, are immediate to calculate explic- 
itly from equations (|4.13p and (|4.14p . We can therefore calculate explicitly all other eigenelements 
by action of Oj^ on those corresponding to the E^,e £ Sq. We summarize this in the following. 

Proposition 5.3. For 7 = (e, ei, 62, • • • e„) e n„>i, set — Ue^Uez ■ ■ ■ Ue„ and Uj = u^^ ° u^.^ o 

• • • . For any 7 G n„>i, we have E^ ~ UjE^i , and 

n 

Ay = Uy(A,0 - A'; A,. + Ai-i . (5.23) 

5.2 Bounded case 

We consider here the case s > d+2. We show that As is bounded and characterize the boundary 
of its spectrum. 

Proposition 5.4. For s > d + 2, Ag is a bounded, and we have 

II A II 1 

\\As\\B{LHdB.dt,)) < C ,d+2-s)/d ' 

i App 

with c = maxjmaxggfg |Ae|,maxeg£ |/3e|}- 

Proof. By Proposition [531 equation (|5.23p . we see that for 7 G n„ we have. 

n 

|A^| < max{max|Ae|,max 

From equation (|5.22p we have As = Ap^*^^-*^'* < 1, therefore the above geometric sum converges, 
and is bounded for all n by its sum 1/(1 — As). □ 
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Wc define the tj-spectrum of as 



Sp^(A,)- fl Sp(A«)\Sp(A,|n„). 

In our case here, this is the boundary of the (pure point) spectrum of A^. Under some conditions 
on A and the /3e, e G £, and for s > d+2 large enough, one can show that Spj^(A,5) is homeomorphic 
to S, and that this homeomorphism is Holder [24]. 

5.3 Weyl asymptotics 

The following theorem justifies calling A^ a Laplace-Beltrami operator. Indeed, for s = sq = d, 
Theorem 15.51 shows that the number of eigenvalues of A^ of modulus less that A behaves like A''/^ 
when A — * cx3, which is the classical Weyl asymptotics for the Laplacian on a compact c?-manifold. 

Theorem 5.5. Let Ms {X) = Card{A' eigenvalue of As : |A'| < A}. For s < d+2, we have the 
following Weyl asymptotics 

C_Xd/(d-s+2) < < ^^^d/{d-s+2) ^ ^5 24) 

as A +00, for some constants < c_ < €+. 

Proof. By Proposition 15. 3i equation (|5.23p . there exists constants x±,y± > 0, such that for all 
7 G n„ we have a;_A^ + < |A^| < x+A"^ + y+. For s < d + 2, A,, > 1 (see equation (|5.22p ). 
so there is an integer k > Q (independent of n) such that a;+A" + yj^ < a:;_A"+'^ + y-- Hence 
for all 7 e 11/, ^ < n, we have \X^\ < x+A" + and for all 7 e 11/,? > n + k, we have 
|A^| > x-A""'"*^ + y_. Therefore we get the inequalities: Card n„ < J\f{x+A^ + y+) < Card n„+i;. 
There are constants ci > C2 > such that for all I G N, ciApp < Card 11/ < C2App, so that 
we get ciApp < 7V(a;+A" + y+) < C2A'^p = csApp. We substitute A^ from equation (|5.22p to 
complete the proof. □ 



5.4 Seeley equivalent 

For the case s = sq = d wc give an equivalent to the trace Tr (e*'^''), for s = sq = d, as t I 0. 
The behavior of Tr (e*'^'') like t^"^/^ as t ], is in accordance with the leading term of the classical 
Seeley expansion for the heat kernel on a compact d-manifold. 

Theorem 5.6. There exists constants c+ > c_ > 0, such that as t I 

^_^~d/2 < rj.^ (gtA,^ < c+t-'^^^ . (5.25) 

Proof. Let P^ be the spectral projection (onto E~f) for 7 G n„>i, and Pq that on Eq. The trace 
reads 

oc 

Tr (e*'^') = 1 + e^''*Tr (Po) + ^^-^^ ' (^'^S) 

n=l 76n„ 

Now Tr (P-y) ~ n-y — 1, with the number of possible extensions of 7 one generation further 
(see equation (|4.14p in theorem 14. 3p . and Tr(Po) = n-o — 1. Since the Bratteli diagram of the 
substitution is stationary, the integers are bounded, so there are p^,p^ > 0, such that for all 
7 G n we have: 

p_ <Tr(P^) <p+. (5.27) 

By equation (|5.23p . the X^Pe, being bounded, there exists A_, A+ > 0, such that for all 7 G 11 we 
have: 

X-Alf' <\X-,\< X+Al%^ . (5.28) 
The cardinality of n„ grows like A'^p so there arc 7r_, 7r_|- > 0, such that for all n > we have: 

TT-K-'pp <\an\<T^+A''pp. (5.29) 
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We substitute inequalities ([07|) . ((05|) . and ((05)) into equation ([g?^ to get 

OO OO 

1 + p-TT^ J2 Apj^e-*^+<^^^' < Tr (e*^=) < 1 + p+n+ ^ A'^^e^*^-<^^' . (5.30) 

n=0 n=0 

Set Nt = (ilog(l/i)/(2 log(Api?)), and split the above sums into two parts: the sum over n < Nt, 
and the remainder. For the finite sum we have: 

Nt-l Nt-1 

J2 A?>pe-*^±<-'' = t-^/2A-^* A^pe-^±<"-^"'""' , (5.31) 

where we have used Ap^^^ = 1/t. With the inequalities e~^^ < e^^±^J=F '' < 1, the above 
sum on the right hand side of (|5.3ip is bounded by the geometric series X)!^=o^ ^^pf ~ i-^PF ~ 
l)/(Api? - 1). Multiplying by t^'^l'^k"^^ we get the inequalities: 

7Vt-l 

i-d/2^/_ < ^ A^pe-*^±^p'-'' < t-'l^d^ , (5.32) 

ri=0 

for some constants c'_, c'_^ > 0, and t small enough. 
For the remainder of the sums in (|5.30p we have 

OO OO 

^ AJpe-*^*^-'-'' = A% A'j?pe-^±Ap"'' = r'^'c'i , (5.33) 

n—Nt m— 

where c'^ > is the sum of the absolutely convergent series. We put together inequalities (|5.32p 
and equation (|5.33p into inequalities (|5.30p to complete the proof. □ 



5.5 Eigenvalues distribution 

We now restrict to the case s = sq = rf, so that A^ = ^pp equation (|5.22p . We also suppose 
here that the weights are simply given by the d-th roots of the measures: diam[7] = /i[7]^/''. 

Let us consider the cases d = \ and d = 2 first. One sees from equations (j4.12bp . (|4.13p . and 
(|5.20p . that the X^,e G and /3e,e G £, are rational functions of App. Hence by equation (|5.23p 
all the eigenvalues of A^; belong to the field Q[App]. The maps in equation (|5.22p become afhne 
maps in this field, whose linear part is the multiplication by Ap''^ (= App or App). 

Let Pa(X) = AT'" + Or-iAT''^^ + • • ■ aiX + ao be the minimal polynomial of A. The field 
Q[App] is isomorphic to , so that the eigenvalues of can be represented by points in Q^. 
The Ae,/3e, are identified with fixed vectors X^,Pe G Q'^, and the maps Ue become afhne maps 
in Q''. The multiplication by Ap''p (linear parts of the maps u^) is implemented in the basis 
Q © Ap/pQ © • • • A%~^^^'^Q by the companion matrix of Pa: 

-ao 

1 ■•. : 

-~ar-2 

1 -a,._i 

Note that it is equivalent to A and therefore has the same Perron-Frobenius eigenvector. Since 
there are finitely many Pe, G : upon multiplication by a large enough integer, we can actually 
represent the eigenvalues by points of the lattice Z''. It is then natural to ask which points 
correspond to eigenvalues of Ag and if one can characterize this set. The answer is surprisingly 
simple: 
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Eigenvalues correspond to points of integer coordinates in Q*" that stay within a bounded strip to 
the Perron-Frobenius eigenline of A. 

See section [01 for an example, and Figure [3] in section [TJ 

For d ^ 1,2 a similar result hold, but the points can no longer be chosen to have integer 
coordinates. In order to prove this in general - for d generic - we consider a quotient ring of 
Q[Api?] where we can implement the multiplication by ^pp- Define Qa{x) ~ Pa{x'^^'^) if d is 
even, and let C be its companion matrix. If d is odd define Qa{x) = Pa{x'''), and let C be the 
square of its companion matrix. Let us denote by d' either d/2 if d is even, or d if it is odd. The 
matrix C implements the multiplication by Ap^ on the ring 'Q\X]/Qa, which we embed in the 
vector space ^ = R © KppM. ® K],pR • • • K%~^R. 

Because of the factor diam[77]^~^ = /i [77]*^^ in equation (I4.12bp . we have to work over 
the field R. An eigenvalue of A^; can then be written as a vector e V . The action of the 
Cuntz-Krieger algebra on such an eigenvalue, given in equation (|5.22p . becomes here 

Me(A^) = CA^ + /3e, (5.34) 

where /3e is the expression of [3^ as a vector in V . The general expression of an eigenvalue given 
in proposition [231 equation (|5.23p . takes here the form 

n 

A^ = C"Ae' + ■ (5-35) 

We now assume that C is Pisot, i.e. it has a single real eigenvalue /i„ > 1, and all its other 
eigenvalues fi satisfy \fi\ < 1. Let us denote by Vu the eigenspace of fiu, and by Vg its orthogonal 
complement in V. We can now characterize the distribution of the eigenvalues of A^. 

Theorem 5.7. There exist a constant m > such that for any eigenvalue X-y of Ad, one has 

distv(A^,K) < Y^^, (5.36) 

where Cs is the projection of C to Vs, and disty the induced Euclidean distance on V. 
Proof. The distance disty (A^, Vu) equals the norm of the projection of A^ to Vs'- 

distv(A^,K) = \\P-^{CsP-^K' +T.%iCi-^P~^Pe,)\\ 

< \\P-'r max{||;3e||,||Ae|| : eefo,ee£} T.%o\\Ci\\ 

< ™ ET=^ liaip- = m/(l-||C.||), 

where P is an invertible matrix diagonalizing C, m ^ max{||/3e||, ||Ae||}, and the series 

converges because Cg is diagonal with eigenvalues of moduli strictly less than one. □ 

Remark 5.8. IfC is no longer Pisot &Mt strictly hyperbolic (i.e. has no eigenvalue of modulus 1), 
then the above result still holds but with Vu the unstable space ( the span of the eigenvectors with 
eigenvalues > 1) and Vs its orthogonal complement. 

6 Examples 

We illustrate here the results of sections |3] and [5] for the classic examples of the Thue-Morse, 
Fibonacci, Ammann-A2, and Penrose tilings. 
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6.1 The Fibonacci diagram 



The Bratteli diagram for the (uncoUared) Fibonacci substitution 



ab 



reads 




where the top vertices are of "type b" and the bottom ones of "type a" , where a = l/cj) — ( VS— 1) /2 
is the inverse of the golden mean, and the term a" at a vertex is the measure of the cyhnder of 
infinite paths through that vertex. 

Note that this substitution does not force the border, so that dB is not the transversal of 
the Fibonacci tiling space. For illustration purposes it is however worth carrying this example in 
details. We treat the "real" Fibonacci tiling together with the Penrose tiling in section [6731 

Since the Bratteli diagram has only simple edges, as noted in Remark 12. 6[ the paths can be 
indexed by the vertices they go through. The paths in IIi are thus written a and &, and the paths 
in 112 are written aa, ab, and ba (note that these are not orthonormal bases for the dot product 
given by the Dixmier trace , so that the Laplacians written below will not be symmetric). The 
restrictions of the Laplace operator (|4.12|) for s = sq = 1 to Hi and 112 ai'c given below together 
with their eigenelements. 



A|n, - 



1 



-1 



with eigenelements (0 , 



) , (1 + 



1 



-I + 302 

-1 + 02 



-1 + 202 
2-402 

1 



-1 + 02 
-1 + 02 



with eigenelements 



(0, 



), (1 + 202 



1 
1 

1-02 



), (3 + 1 



1 

1-02 





)■ 



Using the identities Xb — Xfca, and Xa ~ Xaa + Xa6, we see that the first two eigenvectors of A|n2 
are exactly those of A|ni expressed in 112. Note that 1 — 02 = —0 and that the above eigenvectors 
of A are Xa + Xt, Xa - (f>Xb, and Xaa - (f>Xab- And all other eigenvectors are given by x-yaa - (f>X-(ab 
for 7 e n. 

Since the Bratteli diagram has only simple edges, as noted in Remark 15.11 we can take B ~ B 



and A = A = 



1 1 
1 



The action of the two Cuntz-Krieger operators Ua and Ub on the 
eigenvalues of A as in equation (|5.19p is given in here by 

Ua{X-y) ~ 02 + 1 — 02 , Ub{Xj) = (p^Xy — 1+02 

if 7 is compatible with their action, and Ua(A^) = or Ufc(A^) 

Over the ring Z © 02Z, the companion matrix of A is 
become the affine maps 



= otherwise 

1 

1 1 



(6.37) 

, and the operators (|6.37p 



1 

1 1 



U{,(A, 



1 

1 1 



A. 



-1 
1 
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when 7 is compatible with the corresponding action. 

Figure [3] iUustrates Theorem 1 5 . 71 that characterizes the repartition of the eigenvalues of — A as 
point of integer coordinates that stay within a bounded strip to the Perron- Froebenius eigenline 
of A (slope 0) in Z © cfPZ. Note that the repartition of points in the strip is not "homogeneous" , 
i.e. the number of points within a distance r to the origin is not linear in but rather follows the 
Weyl asymptotics in ^/r (Theorem I5.5p . 



6.2 The dyadic Cantor set and the Thue— Morse tiling 

Those examples have enough symmetries to allow easy and direct calculations (without using the 
operators of the Cuntz-Krieger algebra) . The Bratteli diagram B of the dyadic Cantor set is the 
dyadic odometer, 



B 



and its associated diagram for its Cuntz-Krieger algebra B is the Bratteli diagram of the (uncol- 

^ 01 



fared) Thue-Morse substitution 



1 



10 



B 




where the term ^ at a vertex is the measure of the cylinder of infinite paths through that vertex. 
The top vertices are of "type 0" , the bottom ones of "type 1" . 

We label the paths in B by sequences of O's and I's labeling the edges they go through from 
the root: 7 G n„ is written 7 = (ei, • • • ,e„). The Laplacian on B commutes with the following 
operators: 

^{S::::::r-" S^if"'"' 

where the addition is taken mod 2. The operators Ti commute with each other and square up to 
the identity. One can therefore choose an eigenbasis for A made of eigenelements of the 's: that 
is Haar functions on the dyadic Cantor set dB. We recover this way the example treated in [36] 
and we refer the reader there for the details. 

For the Thue-Morse diagram B, we can also index paths by sequences of O's and I's labeling 
the vertices they go through from the root. The Laplacian is also commuting with the operators 
fi defined like the defined in equation (|6.38p . A basis of eigenvectors of A for s = sq = 1 is 
given by the constant function Xgs ('^ith eigenvalue 0), and the functions 

for n G N, with eigenvalues A„ = — 1(7-4"^^ — l) of degeneracy 2"~^ = Cardn„. The eigenvalues 
satisfy the induction formula A„+i = 4A„ — 2. 



The Weyl asymptotics of Theorem [53] reads here ^J^X + ^ < JV{X) < Jf A + f 
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6.3 The Penrose tiling 

The Fibonacci, Penrose, and Ammann-A2 [21] tilings have formally the "same" substitution on 
prototiles modulo their symmetry groups, with Abelianization matrix 



A = 



2 1 
1 1 



The Penrose and Ammann-A2 substitutions force the border. And for the Fibonacci tiling, one 
considers the conjugate substitution a — > baa, b — > 6a, which is primitive, recognizable, and forces 
the border as noted in Example 12.171 Those three substitution tilings have the same Perron- 
Frobenius eigenvalue, namely Kpp = 0^, where (j) = (l + '\/5)/2 is the golden mean. In conclusion, 
the transversals of those tiling spaces can be described by the set of infinite paths in the same 
Bratteli diagram B illustrated below: 



B 




where a ~ 1/ (f) ~ {sf^ — l)/2 is the inverse of the golden mean, and the term a"' /\G\ at a vertex 
is the measure of the set of infinite paths through that vertex. And where G is the symmetry 
group of the tiling introduced in section [2.241 and \G\ the cardinality of G. That is G = {1} is 
the trivial group (so \G\ = 1) for the Fibonacci tiling, G = C2 y- C2 for the Ammann-A2 tiling 
(symmetries of "vertical and horizontal" refiections, \G\ = 4), and G = Diq for the Penrose tiling 
(10-fold rotational symmetries, and reflections, \G\ = 20). 

Let us denote hy Ui,i = 1, ■■■ 5, the generators of the Cuntz-Krieger algebra (|5.15p associated 
with the Abelianization matrix of B. The induced action on the eigenvalues of A = A^^, as in 
equation (|5.19[) reads here for Penrose and Ammann-A2 (for Fibonacci, A^Jp = 0^ has to replace 
(f)^ in the following equations): 



Ul 



(A,) ^0 A, + (1-0) ^^^^ + ^^^^^ 



.,(A,) = A, + (1 - )^^ + ^^-^ 

.3(A,) = 0^A, + + ^^4(^ + ^^^^ffe^ 

.3(A,)^0A, + (l-0) ^^^^ +^(^ 

if A^ is compatible with the operators. Here G = G^q as in equation (|4.12b[) . so we have G(o) = 
\G\{\G\-l){{w\y + &) + \G?^^, Gisa) = 2((^)2 + (^)2)+4^^ and G(e,) = 2^^. 
The eigenelements of Ajn^ arc for XdB, Aq for Xe„ - ^Xe^: A^^ for Xe„e - <^Xe„/ , e, / G cxti(£a), 
and Ajj, for Xe^e ^ : e, / G cxti(e,,), where 



Ao 



-2|G|(|G| + l-40^) 
|G|2-10|G| + 5 ' 



G(c 



Gisa)' G(o) GiSb) 
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